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FOREWORD 


This  is  the  second  in  a  series  of  reports  on  high-speed 
tensile-impact  work  performed  in  the  Pioneering  Research 
Laboratory  of  NLABS.  The  first  report  (NLABS  Tech.  Rept 
69-18-PR)  appeared  over  5  years  ago.  It  had  been  planned  to 
publish  all  three  of  the  projected  reports  at  that  time,  but  for 
various  reasons  the  second  and  third  reports  were  delayed. 

The  first  report  dealt  with  the  basic  mathematical 
techniques  and  their  application  to  a  linear  material.  The 
present  report  deals  with  non-linear  materials  which  do  not 
exhibit  creep  or  other  time  dependence.  The  third  report  is 
planned  to  deal  with  non-linear  materials  that  do  exhibit  creep. 
Since  all  actual  materials  exhibit  creep  it  is  only  in  the  third 
report  that  we  can  hope  to  demonstrate  good  agreement  between 
theory  and  experiment. 

Body  armor  must  be  resistant  to  high-speed  impact  in 
order  to  stop  projectiles  and  shell  fragments.  Parachute  straps 
are  subjected  to  sudden  tensile  strains  when  the  parachute  opens 
and  the  foiling  load  is  suddenly  decelerated.  The  mathematical 
techniques  described  in  the  present  paper  should  be  useful  in 
finding  out  more  about  how  materials  fail  under  rapid  impact. 
With  this  knowledge,  better  materials  can  almost  certainly  be 
developed; 
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ABSTRACT 


An  integral- equation,  successive- substitution  solution 
for  the  propagation  of  strain  waves  in  a  linear  elastic  material, 
developed  in  an  earlier  paper,  is  here  extended  to  non-linear 
materials  without  creep  (time- independent  materials).  A 
previously- published  problem,  based  on  experimental  data  for 
nylon  string  is  solved  by  the  successive-substitution  method  and 
compared  with  the  earlier  solution  obtained  by  the  method  of 
characteristics?  The  agreement  of  the  two  solutions  is  at  first 
good,  but  in  the  region  where  the  stress-strain  curve  of  the 
nylon  is  concave  upward  a  tendency  to  oscillation  builds  up  and 
the  solution  eventually  oscillates.  This  occurs  in  approximately 
the  same  region  where  shock  waves  can  appear  in  a  method-of- 
characteristics  solution.  Oscillations  in  the  strained  portion  of 
a  string  are  examined  theoretically  and  it  is  shown  that  the  energy 
of  standing  waves  can  be  used  to  obtain  energy  conservation. 
Although  the  successive- substitution  method  generates  oscillations 
when  it  has  to  deal  with  a  discontinuity  in  strain,  the  possibility 
that  oscillations  also  actually  exist  as  a  means  of  conserving 
energy  cannot  be  ruled  out. 
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1.  Introduction 


The  present  report  is  the  second  in  a  series  of 
three  that  describe  the  results  of  mathematical  and  experimental 
investigations  of  the  response  of  polymers  to  tensile  impact* 

Professor  Crout  is  responsible  for  the  mathematical  develop¬ 
ments  and  Mr.  Pilsworth  is  responsible  for  the  experimental  data* 

The  situation  when  the  work  was  undertaken  and  the  reasons  for 
adopting  the  present  approach  are  discussed  at  length  in  the  first 
report  [l]  (subsequently  referred  to  as  Report  I)  of  this  series 
and  will  be  only  briefly  summarized  here. 

In  Report  I  the  mathematical  method  was  developed  and 
was  applied  to  the  problem  of  strain  propagation  in  a  linear 
material  having  no  creep  or  other  time  dependence.  Since  the 
solution  of  this  problem  is  known,  and  since  the  new  method  gives 
the  correct  solution,  its  validity  for  this  case  has  been  established, 
and  by  inference  the  new  method  can  also  be  expected  to  give 
correct  answers  to  related  problems  for  which  satisfactory  solutions 
have  not  been  available. 

The  present  report  extends  the  mathematical  analysis  to 
include  non-linear,  time- independent  materials  and  in  the  third 
report  non-linear,  time- dependent  materials  will  be  treated.  The 
subject  of  the  present  report  (non-linear,  time- independent  materials) 
has  been  treated  by  the  method  of  characteristics  by  Von  Karman  and 
Duwez  [2],  by  G.  I.  Taylor  [3],  and  by  Rachmatulin  [4],  Hence  the 
mathematical  results  of  the  present  report  are  not  presented  as  a 
new  solution  but  rather  as  a  new  method  of  obtaining  a  known  solution 
and  as  the  next  logical  step  toward  solving  the  problem  of  strain 
propagation  in  a  non-linear,  time- dependent  material.  Since  the 
present  method  furnishes  an  alternative  to  the  method  of  character¬ 
istics,  a  comparison  of  the  results  obtained  in  solving  a  particular 
problem  by  the  two  methods  is  given. 

2.  Summary  of  Report  I 

A  satisfactory  theory  of  the  response  of  materials  to  loading 
(tensile  impact)  should  permit  us  to  calculate  the  stresses  and  strains 
in  a  material,  for  any  rate  of  loading,  from  a  few  fixed  properties  or 
parameters  of  the  material.  There  is  no  satisfactory  theory  that  will 
do  this.  Linear  viscoelastic  theory  is  adequate  for  small  strains  but 
not  for  large  strains.  Stress-strain  curves  plus  creep  or  relaxation 
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data  are  satisfactory  for  large  or  small  strains,  if  the  strain 
rate  is  very  low.  The  method  of  characteristics  is  adequate  for 
non-linear  time- independent  materials,  but  as  has  been  shown 
by  Pilsworth  and  Hoge  [5],  by  Pilsworth  [6],  and  by  Smith  and 
Fenstermaker  [7]  it  is  inadequate  for  materials  such  as  Nylon 
and  rubber. 

In  Report  I  a  new  method  of  handling  the  problem  of 
strain-wave  propagation  was  developed,  in  which  the  relevant 
partial  differential  equations  are  transformed  into  integral 
equations.  The  integral  equations  are  then  solved  by  the  method 
of  successive  substitutions.  Referring  to  the  system  shown  in 
Fig,  1,  the  differential  equation 

^  _  2 
a*1  “  **  .ax1 

was  derived  strain  propagation.  Here  x  is  the  position  of 

a  particle  of  material  in  the  unstrained  state,  s  is  the  displace¬ 
ment  of  the  same  particle  from  its  original  position,  and  c  is  the 
wave  velocity.  Other  quantities  are  defined  in  the  table  of 
nomenclature.  By  integrating  eq  ( 1)  twice,  the  equation 

-d.  r-p(x)  4  t  3OO  +  W  j  cttr  (2) 

•'n  "'/1 


was  obtained;  this  equation  is  valid  when  c  is  constant  ( a  linear 
material  has  c  constant).  The  function  f  (x)  specifies  the  initial 
displacements  and  g(x)  specifies  the  initial  velocities;  except 
at  the  impacted  end  of  the  string  these  functions  are  usually  zero 
in  practical  problems. 

Equation  (2)  was  simplified  by  putting  f  and  g  =  o  ahd 
introducing  three  new  normalized  variables: 

-  a*  x  *£  w  2^1 
=  T^o  5  *  "  ~C  ’  *  =  L  (3) 


in  which  v  is  the  velocity  of  the  impacted  end  of  the  string, 
o 


2 


3 


The  result  of  these  substitutions  (see  Report  I,  p  36)  was  to 
reduce  eq  (2)  to  the  form 


-  ff 


cbt  dfc 


where  for  convenience  the  bars  have  been  dropped  from  the  new 
variables.  The  velocity  c  does  not  appear  in  eq  (4) ;  this  equation 
refers,  in  fact,  to  the  problem  in  which  L  =  1,  c  =  1,  and  v  =  1, 
When  a  solution  of  eq  ( 4)  has  been  obtained,  the  results  can  be 
applied  to  actual  problems  by  using  eq  (3),  The  method  devised 
for  evaluating  this  integral  was  to  accept  a  lattice  of  points  in  the 
(x,  t)  plane,  to  specify  the  values  of  s  along  three  borders  of  this 
lattice,  and  to  find  the  values  of  s  at  the  remaining  points  by  an 
iterative  procedure  based  on  eq  (4),  To  find  values  of  the  inte¬ 
grand  the  three- point  differentiation  formula  was  used,  and  for  the 
integration  the  trapezoidal  formula  was  adopted.  Values  of  s  were 
calculated  one  row  (one  value  of  t)  at  a  time  by  iteration,  using  the 
values  of  the  previous  row  as  starting  values.  When  the  specified 
degree  of  convergence  had  been  reached,  the  next  row  of  values 
was  calculated. 


It  was  demonstrated  that  convergence  was  much 
accelerated  by  modifying  the  successive  substitution  process.  The 
modification  adopted  was  to  multiply  the  output  of  the  k  th  iteration 
by  0,  7  and  the  input  by  0,  3,  and  add  the  two  products  to  obtain  the 
input  for  the  next  iteration.  The  weighting  factors  0,  7  and  0.  3 
appeared  to  be  near  the  optimum  for  general  use  but  do  not 
necessarily  give  the  most  rapid  convergence  under  all  circumstances. 
It  was  shown  that  the  use  of  weighting  factors  affects  only  the  rate  of 
convergence,  not  the  final  result. 

Two  linear,  time-independent  problems  were  solved 
by  digital  computer,  using  the  methods  summarized  above,  and 
the  results  have  been  given  in  Report  I.  In  the  first  problem  the 
end  of  the  string  is  set  in  motionwith  its  full  velocity  vQ  at  t  =  o; 
in  the  second  the  velocity  rises  linearly  during  a  finite  time  interval 
beginning  at  t  =  o.  Since  the  exact  solution  of  eq  (4)  is  known,  the 
results  obtained  by  the  iterative  process  could  be  checked  by  compari¬ 
son.  For  the  case  of  a  discontinuous  rise  in  velocity  (infinite 
acceleration  at  t  =  o)  the  values  of  s  found  by  iteration  were  found  to 
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agree  quite  well  with  the  exact  solutions,  but  when  values  of  the 
strain  e  (=  Ss/S  x)  from  the  iterative  solution  were  compared 
with  the  exactTvaiues  it  was  found  that  the  iterative  values  oscillated 
above  and  below  the  exact  values,  when  plotted  versus  x.  The 
oscillations  were  also  present  in  the  values  of  s,  of  course,  but 
were  very  much  less  noticeable  in  s  than  in  9  s/ 3x.  The  oscillations 
were  largest  at  the  values  of  x  nearest  to  the  discontinuity  of  the 
exact  solution.  It  is  clear  that  the  iterative  process  cannot  give 
results  that  reproduce  a  discontinuity  exactly.  This  had  been 
anticipated.  In  adjusting  to  the  discontinuity  the  iterative  results 
were  found  to  be  first  too  low,  then  too  high,  and  so  on,  but  as  the 
discontinuity  (wave  front)  traveled  away  from  any  particular  point, 
the  amplitude  of  the  oscillations  near  that  point  was  found  to  decrease. 

The  case  of  a  continuous  rise  in  velocity  (finite 
acceleration)  was  treated  because  it  corresponds  more  nearly  to 
the  situation  in  actual  experiments,  and  because  it  was  expected 
that  the  iterative  solution  could  follow  the  exact  solution  in  this  case 
with  less  tendency  to  be  thrown  into  oscillation  than  in  the  case  of 
infinite  acceleration.  The  results  justified  this  expectation.  Figure 
2,  reproduced  from  Report  I,  shows  the  strains  calculated  from  the 
exact  solution  (solid  line)  and  from  the  iterative  solution  (dots) 
when  the  velocity  of  the  impacted  end  of  the  string  rises  linearly  from 
0  to  v  as  t  increases  from  0  to  0,  5,  ( Since  the  t  used  in  the  calcula- 

p 

tions  is  the  normalized  variable,  t  =  0,  5  is  the  time  at  which  the  lead¬ 
ing  edge  of  the  strain  wave  has  reached  the  middle  of  the  string, ) 

The  oscillations  of  the  solution  obtained  by  the  successive-substitution 
method  are  noticeable  but  not  large. 

It  is  important  that  oscillations  caused  by  discontinuities 
should  not  destroy  the  validity  of  the  successive- substitution  method. 
Hence  the  effect  of  other  factors  on  the  oscillations  due  to  discon¬ 
tinuities  was  examined.  Since  the  eventual  object  of  the  present  work 
is  to  be  able  to  calculate  the  behavior  of  a  non-linear,  time- dependent 
material,  the  effect  of  non-linearity  was  considered.  When  the  stress- 
strain  curve  is  concave  downward,  the  slope  of  a  strain  wave 
decreases  in  steepness,  and  hence  the  tendency  toward  oscillations  is 
reduced;  but  when  the  curve  is  concave  upward  the  strain  wave 
increases  in  steepness  so  that  there  is  a  possibility  that  a  shock  wave 
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Strain  for  a  linear  time-independent  model  as  obtained  from  digital  computer  calculations 
ith  j  =  105  (t  =  1.  75).  The  velocity  of  the  impacted  end  increased  linearly  to  v  at  t  =  C 


will  form.  Since  a  shock  wave  is  a  discontinuity  it  will  increase 
the  tendency  toward  oscillations  in  the  iterative  solution.  Time- 
dependence  of  the  material  will  tend  to  delay  a  part  of  the  response 
and  hence  will  tend  to  destroy  discontinuities  and  reduce  the  tendency 
toward  oscillations.  The  conclusion  reached  in  Report  I  was  that 
the  modified  successive  substitution  procedure  was  satisfactory;  and 
that  the  oscillations  relative  to  the  exact  solution  could  be  trouble¬ 
some  but  would  probably  not  be  serious,  except  possibly  where  the 
stress- strain  curve  was  concave  upward  and  shock  waves  were 
mathematically  predicted. 

3.  Extension  of  the  Solution  to  the  Non-Linear  Case. 

The  strain  at  any  point  in  the  string  is  e  =  0s/ fix. 

In  the  non-linear  case 


<j~  -  (T  (6  )  j  (5) 


and  Newton* s  law  yields  the  differential  equation 


^  £ 
f  2* 

(6) 

The  integral  equation  is 

ft  t 

t)  -  £  get)  *  j J J  -|T  m dt 

(?) 

where 

al(x,o)  -f(y-)-°  ~) 

^  &U-.0  J 

(8) 

Since 

(9) 

and  since  the  wave  velocity  c  is 
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it  follows  that  ( 7 )  may  be  written 


Jt  £ 

t)~£  g(x)+£j f  (“) 

It  should  be  noted  that  c  is  now  a  function  of  £  ,  which  is 
8  s/  gx,  Let  us  place 


where 

0 

which  is  the  wave  velocity  for  zero  strain* 
Substituting  in  (11)  we  obtain 


or 


J,  (.%,£) -i 


(15) 


In  the  present  case 

^  (X  )  =  O  0^  *  ^  I  ) 


(16) 


it  being  noted  that  this  relation  applies  at  the  right  end  ( no 
discontinuity  in  velocity) . 
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Equation  (  15)  thus  reduces  to 


Ci,*) 

c  o 


(17) 


Here  c  is  a  function  of  £  ;  however,  it  follows  from  (12)  that 


€  =• 


-  ^  da  - 


L 


*  ^  =  e  • 


hence 


(18) 


*.(*■)**.<£)  ^  (|t)- 


(19) 


It  is  now  evident  that  the  technique  involved  in  the  present  report 
differs  from  that  in  Report  I  only  in  that  the  factor  (c/cq)2  must 
be  included  in  the  integrand. 

Finally  we  note  that  in  addition  to  ( 16)  we  have 

Z(o,±)  =  O  o<  f 

Z(x ,0)-  O  oix-i/ 

0-  I^i  j 

also  in  regard  to  velocity 

/V  -  -  —  TZ  sCaiW. 

4.  Description  of  the  Repetitive  Process  When  Applied  to  the 
Non-Linear  String. 

The  following  set  of  instructions  for  programming  the 
successive  repetition  process  for  use  in  connection  with  a  digital 


(20) 


(21) 
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computer  was  devised.  A  new  quantity 


*  ’  ■  0  •  -pg 


(22) 


not  used  in  Report  I  is  needed.  The  weight  factors  a 
as  described  in  Report  I,  p  51,  are  again  used. 


=  .3,  p=  ,7, 


Our  problem  is  to  determine  the  values  of  the  quantity 
s  - .  at  the  points  of  the  square  lattice  shown  in  Fig.  3.  The  index  i 

J 

takes  the  values  0,  1,  2,  3,  ...»  n;  and  j  takes  the  values  0,  1,  2, 

3,  . . .  .  Since  n  is  given,  the  range  of  values  taken  by  i  is  bounded. 

The  range  of  values  taken  by  j,  however,  has  no  upper  bound  except 
one  which  will  be  assigned  later  from  practical  considerations. 


The  values  of  sjj  are  given  for  those  points  which  lie 
on  the  border,  thus 


*  0 

S- 0,1,2, 3  ..  • ; 

-4^0  -  ° 

^  r  6,  1,2,3  •  •  ‘*lj 

(^/oLsyO 

j  =  0,  1,2,1  .  ..  . 

(23) 


The  values  of  s^  at  the  other  points  are  determined  a 
row  at  a  time  —  that  is  for  j  fixed  and  i  =  1,  2,  3,. . . ,  (n-1). 


The  values  of  j  are  taken  in  the  order  1,  2,  3  . . .  .  In  view  of  this 
it  will  evidently  be  sufficient  to  describe  how  to  obtain  the  values 
of  s..  for  the  j**1  row  when  the  values  for  all  of  the  preceding  rows 


ij 


are  known. 


The  values  of  for  the  j111  row  are  obtained  as  the 
limits  gotten,  respectively,  by  applying  repeatedlg^a  process 
shortly  to  be  described.  At  the  beginning  of  the  k  repetition  we 


have  the  quantities 


At 


<*) 


Ok)  Ok) 


'j 


A, 


)  )  3 j  ’  '  '  > 

as  approximations  to  the  quantities 


At 

/U-J,  J 


ik) 


(24) 
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(25) 

;  >  -*aj  >  5  *  '  '  y 

respectively.  Here  k  is  a  superscript,  not  an  exponent.  When 
the  process  is  applied  starting  with  the  quantities  (24),  the 
quantities 


-4>, 


c-fc+o 


(.Jk  +  i)  (Jn-I) 


I  J 


Aj 


•3J 


(■fe  +  l) 


'  J 


■vn-i 


(26) 


are  obtained  as  approximations  to  (25),  respectively.  The 
sequence  of  repetitions  is  started  by  taking  the  values  of  s.. 
for  the  (j-1)  row  as  approximations  to  the  values  for  the 
j  row,  respectively,  thus 


a) 

l  j  m^fj  j-i  > 


( i )  0> 

~  j-i  j  •  •  •  ;  *  ^V-«,  j-i  * 


(27) 


These  values  are,  of  course,  known.  The  sequence  of 
repetitions  is  terminated  when,  for  sufficiently  large  k,  the 
output  (26)  duplicates  the  input  (24)  to  the  required  accuracy. 
The  common  values  which  then  compose  (24)  and  (26)  are 
taken  as  the  final  result  (25).  The  values  of  s*-  which  compose 

x.l_ 

the  jm  row  are  thus  determined. 

Finally  we  shall  describe  the  process  by  which 
we  start  with  (24)  and  obtain  (26).  During  the  course  of  the 
calculations  we  compute  and  record,  in  addition  to  the  quantities 
s--,  the  auxiliary  quantities  £••,  R->,  A^-,  and  Bj..  The  quantity 
A  is  a  measure  of  the  first  time  integral  and  the  quantity  B  of 
the  second  time  integral.  Both  A  and  B  differ  somewhat  from 
the  corresponding  quantities  in  Report  I  because  they  contain  R, 
which  was  not  used  in  Report  L  The  four  auxiliary  quantities 
do  not  exist  at  the  points  of  the  left  border  of  the  lattice,  for 
which  i  =  0,  or  the  right  border,  for  which  i  =  n.  Along  the 
lower  border,  where  j  =  0,  e^,  A.j,  and  By  are  defined  to  be 
zero,  and  Ry  to  be  unity.  As  is  the  case  with  sy  the  auxiliary 
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quantities  £..,  R..,  A..,  and  B..  are  obtained  one  row  at  a  time, 
are  known  in  the  first  (j-1)  rows,  and  are  obtained  in  the  iT1  row 
as  the  limits  approached  by  the  quantities  e/.  ,  R.'.  ,  A.(.  ,  b{.  , 

respectively,  during  the  course  of  the  successive  repetitions,  ^ 
The  desired  basic  process  is  now  described  by  the  relations 

(*)  ^  -fe  \ 

=  SL  /;  Us  -  ‘  •  j  (to*-  I*)  (28) 


R«  -  R  (*fi) , 


=  /,  4,  .  •  ^  (n-l) 


A+f  -  ^,j-l  +  ^4-1  +Ai+Ii4-|)-Ai| i-1  J  ) 

«■*$[* <4$  J  r30) 

(w-  I)  J 

B'*  "  4  J5.{kx,i-<  +•  A1^),  ^ 

,(*♦<)  ,  (to  ,  R.(to 

-  •  3  ■'Va  j  5  •'<-*  -  '  )  A  y  .  . .  j  (/Vt- 1)(  32) 

(k) 

Starting  with  the  given  quantities  s;.  ' ,  j  being  fixed  and  i 
taking  the  values  1,2,  ,,,,  (n-l),\e  compute  the  quantities 
€\.  using  (28);  then  the  quantities  R(k'  using  (29).  The 
function  R  (e)  is  given,  it  is  ^ 


ju&iMj  (ir^srtA^  cOtL  & 

P  (c)  -  - —-3 - ~Z -  ,  (3- 

and  is  expressed  in  a  form  which  is  convenient  for  machine 
calculation—  an  algebraic  polynomial  for  example.  After  R 
has  been  calculated  or  otherwise  found,  the  quantities  A<.k)  ar 
computed  using  (30);  then  the  quantities  B>.  using  (31^j  and 
finally  the  quantities  s(k+1)  using  (32).  T^iis  process  is  now 
repeated,  and  so  on. 


is  now 
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5.  Results  of  Digital  Computer  Calculations  for  a  Non-Linear 
Model 


Experimental  measurements,  made  on  Nylon  by 
Pilsworth,  have  been  previously  published  [5,  6].  Some  of 
the  published  data  had  been  analyzed  by  the  method  of 
characteristics,  and  it  was  decided  that  the  results  of  this 
analysis  would  be  convenient  for  comparison  with  results 
obtained  by  the  present  method  of  successive  substitutions.  In 
the  earlier  investigation,  several  different  stress-strain  curves 
were  used  for  method- of- characteristics  calculations;  the  curve 
that  gave  the  best  results  is  shown  in  Fig.  4.  It  was  decided  to 
use  this  curve  as  the  basis  of  the  comparison. 

Table  1  contains  the  values  read  from  the  curve. 

For  reasons  to  be  discussed  later  it  was  impossible  to  use  the 
iterative  procedure  with  these  stress-strain  values  at  the  highest 
strains,  and  in  order  to  obtain  an  iterative  solution  at  the  highest 
strains  the  upper  part  of  the  stress- strain  curve  was  simply 
replaced  by  a  straight  line.  This  line  is  shown  dashed  in  Fig.  4 
and  the  corresponding  values  are  given  in  parentheses  in  Table  1. 

TABLE  1 

Dynamic  Stress- Strain  Curve  for  Nylon  Yarn 


cr,  gms/den 

c,  cm/ ms 

0 

240 

01 

.  55 

192 

02 

.91 

172 

03 

1.22 

159 

04 

1.  50 

153 

05 

1.76 

150 

06 

2(02 

154 

07 

2.35 

175 

08 

2.75 

220 

09 

3.  44 

314 

10 

4.95 

408 

11 

7.33(6,  84) 

502(408) 

12 

10.  75(8.73) 

596(408) 
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Fig.  4.  Dynamic  stress- strain  curve  for  nylon  yarn. 
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Velocity  cm /milUsec 


D  Lsplacement 


The  problem  that  had  been  previously  solved 
by  the  method  of  characteristics,  and  which  is  now  solved  by 
the  method  of  successive  substitutions,  was  that  of  a  50  cm 
length  of  Nylon  string  with  one  end  fixed  and  the  other  impacted 
by  a  slider  that  accelerated  it  from  rest  to  a  velocity  of  90  m/sec 
in  0.  15  millisec.  The  increase  in  velocity  was  not  linear;  the 
values  used  were  determined  from  the  experimental  data.  For 
the  method  of  characteristics  the  value  of  vQ  during  the  period 
of  acceleration  must  be  known;  and  for  the  iterative  method  the 
values  of  displacement,  s,  at  the  moving  end  of  the  string  are 
required.  The  accepted  values  of  vQ  and  s  are  shown  in  Fig.  5. 
The  values  of  S  were  derived  from  the  previously-used  values  of 
vQ  and  are  consistent  with  them. 

The  rise  time  of  the  velocity,  as  determined 
from  the  experimental  data,  was  about  0.  15  millisec;  this  value 
is  typical  of  the  rise  times  usually  observed.  The  velocity  of 
propagation  of  small  strains,  according  to  Table  1,  is 
240  cm/millisec.  At  this  velocity  the  initial  strain  will  travel 
36  cm,  or  72%  of  the  length  of  the  string,  before  the  moving 
end  reaches  its  full  velocity.' 

Figure  6  shows  the  results  of  the  iterative 
(successive  substitution)  solution  of  the  problem.  The  graph 
shows  strain,  e,  as  a  function  of  position,  x.  Each  plotted  line 
corresponds  to  a  fixed  time,  that  is,  to  a  selected  value  of  the 
time-index  j.  A  square  lattice-mesh  with  n  =  50  in  the 
normalized  variables  was  used.  In  the  actual  variables  this 
corresponds  to  an  interval  in  x  of  L./50  =  1  cm  and  in  t  of 
L./50  c0  =  4.  1667  microsec.  In  the  upper  part  of  Fig.  6  there  is 
evidence  of  oscillations;  these  will  be  discussed  later.  For  the 
present  we  shall  confine  our  attention  to  the  lower  part  of  the 
figure  where  no  oscillations  occur. 

The  previously  published  [5,6]  solution  of 
the  problem  as  obtained  graphically  by  the  method  of  character¬ 
istics  is  shown  in  Fig.  7.  The  abscissa  is  x,  the  ordinate  is  t, 
and  the  straight  lines  are  characteristics;  except  in  the  regions 
where  characteristics  cross  they  correspond  to  lines  of  constant 
strain.  Since  the  coordinates  of  Figs.  6  and  7  are  different  they 
cannot  be  readily  compared.  To  permit  a  comparison  a  single 
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Fig*  6*  Strain  pattern  in  nylon  yarn  calculated  by  successive  substitution  for 
a  time- independent  model.  Time  interval  between  lines  is  0.  0417  millisec* 
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new  graph  has  been  constructed  showing  £  as  a  function  of  x, 
and  containing  typical  results  from  both  solutions.  This  graph 
(Fig.  8)  contains  a  few  of  the  lines  of  Fig.  6  replotted,  and  in 
addition  some  corresponding  values  (plotted  points)  obtained 
from  the  method-of-characteristics  solution.  The  agreement 
between  the  points  and  the  curve  is  quite  satisfactory;  the  dis¬ 
crepancies  are  no  larger  than  might  be  expected  when  one  of 
the  methods  is  graphical.  The  two  points  near  the  top  of  the 
figure  that  do  not  fall  near  a  curve  lie  in  the  region  where  the 
method  of  characteristics  indicated  the  presence  of  shock  waves 
and  the  computer  solution  gave  serious  oscillations.  Since  at  all 
earlier  times  the  iterative  method  gave  the  same  result  as  the 
method  of  characteristics,  its  validity  has  been  established  for 
the  solution  of  a  non-linear,  time- independent  problem,  so  long 
as  oscillations  in  the  solution  are  not  encountered. 

We  now  return  to  a  discussion  of  oscillations 
that  occur  in  the  iterative  method  in  the  region  of  large  strains. 

The  first  evidence  of  oscillations  in  Fig.  6  occurs  at  j  =  110, 
when  the  strain  in  the  material  varies  from  about  0.  051  to  about 
0.094.  Referring  to  the  stress-strain  curve.  Fig.  4,  we  see  that 
it  is  concave  upward  throughout  almost  all  of  this  range  of  strains. 
When  the  stress-strain  curve  is  concave  upward  we  may  (from 
the  mathematical  standpoint  at  least)  expect  shock  waves  to 
develop.  Since  a  shock  wave  approximates  closely  to  a  discontinuity 
we  should  expect  it  to  throw  the  iterative  solution  into  oscillations; 
it  was  shown  in  Report  I  that  oscillations  are  generated  when  the 
iterative  solution  attempts  to  conform  to  a  discontinuity. 

It  may  be  asked  why  the  iterative  solution  does 
not  generate  oscillations  before  j  =  110,  since  at  that  point  almost 
all  of  the  sample  is  within  the  region  in  which  shock  waves  should 
develop.  A  reasonable  explanation  of  this  is  that  the  finite  rise 
time  of  the  velocity  of  the  moving  end  of  the  string  generated  a 
sloping,  not  a  steep,  wave  front  traveling  down  the  string.  Further¬ 
more,  the  initial  portion  of  the  stress- strain  curve  is  concave  down¬ 
ward,  and  this  flattens  the  wave  front  still  more.  As  soon  as 
the  part  of  the  curve  that  is  concave  upward  is  reached,  the  wave 
front  begins  to  steepen  but  only  at  j  =  110  has  it  become  steep 
enough  to  generate  oscillations. 
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Fig.  8.  Strain  pattern  for  nylon  yarn  calculated  by  successive  substitution 
(lines)  and  by  method- of- characteristics  (points). 
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Between  j  =  110  and  j  =  132  the  iterative, 
oscillatory  solution  continued  to  converge,  but  with  the  oscilla¬ 
tions  increasing  in  amplitude.  At  j  =  133  the  oscillations 
rapidly  increased  beyond  the  capacity  of  the  computer.  Modifi¬ 
cations  were  then  made  in  the  program  to  try  to  reduce  the 
oscillations  at  least  to  the  point  where  the  program  could  be 
completed.  Decreasing  the  lattice  spacing  reduced  the  wave¬ 
length  of  the  oscillations  but  didn*t  reduce  their  amplitude 
significantly.  Varying  the  error  criterion,  that  is,  the  maximum 
allowable  difference  between  the  results  of  two  successive  sub¬ 
stitutions  when  the  values  were  accepted  as  converged,  was  no 
help.  Using  a  five-point  least  squares  formula  instead  of  a 
three-point  formula  to  calculate  strain  from  displacement 
reduced  the  oscillations  somewhat  but  only  delayed  the  breakdown 
from  row  133  to  row  137.  Finally,  it  was  decided  to  modify  a 
part  of  the  stress- strain  curve  that  was  being  used,  the  solid  curve 
of  Fig.  4.  A  large  section  of  this  curve  is  concave  upward,  and  a 
part  of  this  section  was  replaced  by  a  straight  line.  The  modified 
part  of  the  stress- strain  curve  is  shown  by  the  dashed  line  in 
Fig.  4.  Examples  of  the  results  computed  from  the  modified 
stress- strain  curve  are  the  curves  for  j  =  140,  150,  and  160  in 
Fig.  6;  oscillations  are  present  but  they  are  relatively  small  and 
do  not  cause  the  solution  to  diverge. 

The  substitution  of  a  straight  line  for  a  portion 
of  the  stress-strain  curve  that  was  concave  upward  was  made 
principally  to  see  what  effect  it  would  have  on  the  oscillations. 

It  is  suggested  by  the  graphical  procedure  that  must  be  used  in 
dealing  with  shock  waves  by  the  method  of  characteristics,  in 
which  the  states  before  and  after  a  shock  lie  at  the  ends  of  a 
chord  joining  two  points  on  the  stress-strain  curve.  Aside  from 
this  the  use  of  the  straight-line  segment  has  no  particular  signifi¬ 
cance.  However,  a  study  of  the  equations  of  motion  in  regions 
where  shocks  were  possible  led  us  to  the  conclusion  that  in  some 
situations  actual  physical  oscillations  might  be  generated,  if  the 
model  we  were  studying  were  sufficiently  sophisticated  to  repre¬ 
sent  an  actual  material.  This  possibility  is  discussed  further  in 
the  next  section. 
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6.  Nature  of  the  Oscillations  when  the  Stress -Strain  Curve  is 
Concave  Upward 


The  oscillations  that  occurred  in  the  solution 
of  the  linear  case  were  explained  in  Report  I  as  computational 
hunting  above  and  below  the  true  solution,  and  they  were  attributed 
to  the  inability  of  the  computational  technique  to  conform  to  dis¬ 
continuities.  So  when  oscillations  are  encountered  in  the  solution 
of  the  present  (nonlinear)  case,  it  is  natural  to  assume  that  they 
are  also  computational.  The  fact  that  the  oscillations  do  not 
become  troublesome  until  the  stress- strain  curve  is  concave 
upward  tends  to  support  this  assumption,  for  when  the  curve  is 
concave  upward  the  front  of  a  strain  wave  becomes  steeper  and 
a  shock  wave  begins  to  develop.  Since  a  shock  wave  is  for  mathe¬ 
matical  purposes  a  discontinuity  our  computational  techniques 
should  have  difficulty  in  conforming  to  it  and  oscillations  should 
be  expected.  We  accept  this  argument  as  valid;  however,  it  does 
not  appear  to  be  the  full  explanation  of  the  oscillations  observed. 

The  much  greater  amplitude  of  the  oscillations 
in  the  non-linear  case  (often  exceeding  the  capacity  of  the  computer) 
led  us  to  consider  whether  the  oscillations  might  not  be  partly 
computational  and  partly  "  real"  .  By  "  real"  we  mean  that  the  exact 
solution  of  the  problem  might  show  oscillations  even  if  the  oscilla¬ 
tions  introduced  by  our  computational  techniques  were  not  present. 

It  should  be  remembered  that  the  problem  whose 
solution  is  shown  in  Fig.  6  is  only  an  approximation  to  an  actual 
problem.  The  problem  solved  in  the  figure  does  not  contain  the 
time- dependency  that  we  believe  is  present  in  actual  polymers. 
Time- dependency  tends  to  smooth  out  discontinuities  and  might 
prevent  oscillations  from  developing  in  an  actual  material.  We 
know  of  no  cases  where  shock  waves  have  actually  been  produced 
by  tensile  impact  on  polymers.  Nevertheless  there  is  some  indica¬ 
tion  that  oscillations  may  be  observed  experimentally.  Much  of  the 
data  taken  in  this  laboratory  appear  to  show  small  oscillations  in 
the  strain  above  and  below  a  smooth  curve  of  strain  versus  position. 
These  oscillations  are  at  the  very  limit  of  the  sensitivity  of  our 
measurements  and  we  have  always  smoothed  them  out  at  an  early 
stage  of  our  analysis. 
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Petterson,  Stewart,  Odell  and  Maheux  [  8] 
also  observed  oscillations  in  strain  patterns.  Their  strains 
were  produced  by  transverse  impact  of  bullets  on  nylon  yarns; 
the  oscillations  were  quite  large.  These  authors  believe  the 
oscillations  were  not  the  result  of  any  defects  in  the  measure¬ 
ment  but  were  actually  present.  However,  since  no  explanation 
of  the  oscillations  was  available,  they  were  smoothed  out  before 
the  data  were  analyzed. 

We  shall  now  present  some  theoretical  con¬ 
siderations  on  the  problem  of  whether  oscillations  should  occur 
in  the  tensile  impact  of  a  time- independent  material  with  a  stress- 
strain  curve  that  is  concave  upward  over  part  of  its  length. 

7.  Theoretical  Investigations  Pertaining  to  Oscillations. 
Momentum  Theorem. 


In  what  follows  we  shall  need  a  momentum 
theorem  which  enables  us  to  investigate  the  momentum  balance  in 
a  region  whose  shape,  size  and  motion  are  specified  arbitrarily. 

At  time  t  let  the  closed  surface  S  in  Fig.  9  contain  a  mass  m  dis¬ 
tributed  in  any  manner,  the  momentum  of  this  mass  being  the  vector 
M,  and  the  vector  sum  of  the  external  forces  acting  on  this  mass 
being  F.  An  infinitesimal  time  dt  later  the  mass  m  has  moved,  and 
occupies  a  volume  bounded  by  the  surface  At  this  time  the 
momentum  has  increased  by  the  amount  dM2  given  by 

dM2  =  F  dt.  ( 34) 

The  volume  whose  momentum  balance  we  are 
investigating  coincides  at  time  t  with  that  bounded  by  the  surface  S. 

At  time  t  it  hence  contains  the  mass  m  and  the  momentum  M.  At 
time  t+dt  this  volume  has  changed  in  a  prescribed  manner  and  is 
bounded  by  the  surface  Sj.  The  momentum  of  the  mass  contained 
within  it  has  increased  by  the  amount  dM1#  which  may  be  written 

dM|  =  dMj  -  dM2+dM2  =  (dMr  dM2)  +  F  dt  (35) 

due  to  (  34) .  But 

I 
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mass.  At  t  +  dt  mass  has  moved  to  S2  and  position  of  S  has  moved  to  S1# 
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dv\,  d  Mx  -  (nrurwurf"^  ^  s,  <cL  Xc«hjl  £+cbb 

—  nrrurn^ctu^  -Uo  S  act  fc) 

—  (  /v>taW/vi£^  as**  act  ctc**&  t  tdt 

~  /^n 4yKOctus»v  js* o  5  HXavJL  t  ) 


-  /m^r>vi/4cAp»^  u>vc  S , 


—  AvuryH^oU^^  -vk. 


cjt  t  +  dt 

^ci*^  £>  i~  dt 


.  S ,  isvdtl  trvrfjuAj*' 

S  2.  oJL  iv+*A*  ts  f'  ^4^  /) 


—  ^  jl^aaJLu  (y  S  ^  ixdd 

^  S/  dt  tri-dt) 

£  o£  />kalaAj  Q/r&MJvuru^  S  ( 

djM^^  ijs»*jL,  Xs*£iw*&  <Jjb  ) 

—  (^/ vri6>M&+tM»+*  a*~A/IAj  *2jux/\tc/y^ 

djXAMi^  (dk  ) 

yyiszt, 

As*-j  ^C4nt-  A +ipWlJL  cfcfc 

JtMJU  zfc  hsJUjxA^y  /^urfZfVV  S/ 

CL^Jl  ~ZU  (36) 
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Substituting  in  (35)  and  dividing  by  dt  we  obtain 

~ Jju*  —  F  JxJxtx,  a/k tALOAUf  ^  AMTM^nJjbiA*-'  1/**  S , 

Amjl -tir  fU  zmJZ^o  *£  s,  (37) 

tS  Azkju  SHU44AJ. 

The  second  term  on  the  right  hand  side  is  of  the  nature  of  a 
convection  term.  Equation  (37)  is  the  desired  theorem.  Using 
it  we  shall  derive  certain  consequences  of  the  assumption  that 
an  abrupt  increase  in  strain  can  be  propagated  at  constant 
velocity  c. 

Propagation  of  an  Abrupt  Increase  in  Strain.  Let  us  suppose  that 
the  strain  increases  abruptly  from  0  to  e  at  a  vertical  wave  front 
which  travels  to  the  left  with  a  velocity  c,  as  indicated  in  Fig.  10. 
As  usual  x  indicates  distance  along  the  string  in  its  unstretched 
state.  Let  us  cbnsider  the  region  bounded  by  the  two  arbitrarily 
chosen,  stationary  dotted  lines,  which  at  time  t  lie  close  to  the 
wave  front,  one  on  each  side,  as  shown.  These  lines  play  the  role 
of  in  (37). 


The  only  external  force  acting  is  cr  A,  where  A 
is  the  cross  section  of  the  (unstretched)  string,  this  being  exerted 
by  the  string  on  the  right. 

In  time  dt  the  wave  progresses  a  distance  cdt  to 
the  left.  The  length  cdt  of  string  thus  undergoes  a  change  in  strain 
from  0  to  e;  hence  a  length  cedt  of  (stretched)  string  passes 
through  the  right  hand  boundary  (dotted  line).  Dividing  by  dt  we 
obtain  for  the  velocity  v  of  the  string  at  this  boundary 


v  =  ce.  (  38) 

Since  the  position  of  this  boundary  was  chosen  arbitrarily  it  follows 
that  the  velocity  is  the  same  and  given  by  (  38)  at  all  points  of  the 
string  to  the  right  of  the  wave  front. 

Let  A  be  the  cross-section  area  of  the  string. 
Then  in  time  dt  a  mass  r  Ac  dt  undergoes  an  abrupt  change  from 
velocity  0  to  velocity  ce.  The  corresponding  momentum  change 
is  c2ePAdt.  It  fdlavs  that  the  increase  in  the  momentum  contained 
in  the  test  region  bounded  by  the  two  dotted  lines  is 
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KA^A &CtZ  Aaa  V j  /ulji  ^  f  A 


C%  6  f 

UxrU*~dAAQ 


). 


(39) 


But  (37)  states  that  this  quantity  is 


F  <£fc  -f 

—  P  ebb  —  AH«Wt«jfcw<.  xA^ui^  ^itc 

^u-  u±U,h*Ak 


(40) 


Equating  (39)  and  (40)  gives 

(<^7'  £  A  ott»  —  F  o(£  ~  <T~  A  djz 

c\e  =  <r 

^-yW  ■ 

Note,  If  the  boundaries  of  the  test  region  had  been  made  to  move 
with  the  string,  the  convection  terms  would  be  missing  from  (39) 
and  (40);  however,  since  these  terms  cancel,  the  result  (41) 
would  remain  unaltered. 


Next  let  us  consider  the  energy  balance  of  the 
test  region  between  the  dotted  lines  shown  in  Fig,  10;  however, 
this  time  we  shall  let  the  boundaries  of  this  region  move  with 
the  string.  In  time  dt  we  have,  noting  (38)  and  (41), 


A  s  >C  €  <T  A  cbt 


(42) 
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(43) 


JMA^(o±**XU*  £**)  ='  pA  JUcUt.(<(f 


-  ~  <.e  <r  A  obt. 


(44) 


Also 


■OHCAjl&AAa  AMj 

(  at  wtML  £ww±)  ~  h  &  cit 


(45) 


where  <r  is  the  average  value  of  cr  (e)  over  the  strain  interval 
ave 

o  to  £. 

The  corresponding  power  equations  are  obtained  by  dividing 
(42),  (44),  and  (45)  by  dt.  We  note  that  exactly  half  of  the 
energy  input  is  utilized  in  increasing  the  kinetic  energy,  which 
leaves  an  amount  z  cecrAdt  to  balance  the  increase  in  potential 
energy. 


In  the  linear  case,  where  cr  is  proportional  to 

e,  it  is  true  that 

<n«  -  -j0' 

and  we  get  a  perfect  energy  balance.  If,  however,  the  stress- 
strain  curve  is  concave  downward 

Avt  ^  ~Z  <r~ 

and  the  available  input  energy  is  not  sufficient  to  provide  the 
required  increase  in  potential  energy.  The  abrupt  wave  front 
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postulated  cannot  be  maintained.  Actually,  the  top  of  the  wave 
lags  behind  the  bottom,  the  wave  spreads  out,  and  the  front  loses 
its  abruptness. 


Finally,  the  case  in  which  we  are  primarily 
interested  is  that  in  which  the  stress-strain  curve  is  concave 
upward,  in  which  case 

sr*  "k  C~~ 

0  avt  Z 


and  the  amount  of  input  energy  available  is  more  than  enough  to 
provide  the  required  increase  in  potential  energy.  The  postulated 
abrupt  wave  front  moving  with  a  constant  velocity  is  hence  not  in 
dynamic  equilibrium,  and  cannot  be  maintained.  In  an  actual 
problem  the  possibility  exists  that  mechanical  energy  may  be 
dissipated  as  heat.  Such  a  mechanism  cannot  explain  the  present 
situation,  however.  The  laws  of  classical  mechanics  must  be 
obeyed,  even  while  kinetic  energy  is  being  converted  into  heat.  It 
may  be  mentioned,  however,  that  the  excess  energy  might  produce 
oscillations  such  as  those  mentioned  on  page  23.  This,  of  course, 
infers  a  modification  of  our  present  abrupt  form  of  wave. 


Effect  of  the  Merging  of  Two  Waves. 

Further  insight  into  the  nature  of  the  difficulties 
encountered  when  the  stress-strain  curve  is  concave  upward  can  be 
obtained  by  considering  two  successive  abrupt  waves,  and  seeing 
what  happens  when  one  overtakes  the  other.  Fig.  11.  Here  the  two 
waves  were  started  at  different  times.  In  time  dt  the  two  waves 
progress  by  amounts  Cjdt  and  Cjdt,  respectively;  the  string  between 
the  two  fronts  moves  a  distance  CjCjdt,  and  the  string  on  the  right 
of  the  second  front  moves  a  distance  [  c^j  +  cz(  €z- £i)  ]  dt.  Dividing 
by  dt  we  obtain  the  corresponding  velocities,  thus 


c*€,  +*/Ca(€t-0=  , 


(46) 


Equation  41  can  be  applied  to  give  the  velocity 
of  the  first  front,  thus 

■c.,  =  Va“i/re'  ri  -  <r<€0- 
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(47) 


In  order  to  obtain  the  velocity  of  the  second  front  we  shall  apply 
the  momentum  theorem  (37)  to  the  test  region  lying  between  the 
two  vertical  dotted  lines  in  Fig.  11,  which  lines  move  with  the 
string.  In  time  dt  the  momentum  in  this  region  increases  by  the 
amount  f»Ac2dt(  V2  -  v^,  the  convection  term  is  zero,  and  the  total 
external  force  is 


f  =  (<r4 -  "1)  A  <r,=  <n&Y 


(48) 


Equation  (37)  thus  becomes 

f  A  (nr2-nr, )-  (<ra  -  <r,)  A 

or  noting  (46) 

•f  cl  Cet-e,)  =  <n  - 


(To  -  <T\ 


fCe»-ei)  -  <49> 

We  are  supposing  that  the  stress- strain  curve 
is  concave  upward;  hence,  we  see  from  (47)  and  (49)  that 


^2  >  /C/  .  ( 50) 

It  follows  that  the  second  wave  is  certain  to 
overtake  the  first.  When  it  does  we  again  have  a  single  abrupt 
wave  front,  the  strain  discontinuity  being  e 2»  If  this  wave 
proceeds*  with  velocity  c,  it  follows  from  (41)  that 


also  the  string  velocity  immediately  to  the  right  of  this  wave  is 


(V  —  /C. 


See  Appendix  2  for  some  comments  on  an  alternative  approach, 
in  which  the  merged  wave  is  assumed  to  propagate  with  a  velocity 
lower  than  c. 
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due  to  (38).  The  difference  quotients  which  appear  under  the 
radicals  in  the  expressions  for  Cj,  Cz»  and  c  are  the  slopes  of 
the  three  chords  shown  in  the  stress-strain  curve  of  Fig.  12. 

Just  before  the  second  wave  catches  up  with 
the  first  the  velocity  of  (almost)  the  entire  string  to  the  right 
of  the  wave  is  V2.  Just  after  the  second  wave  overtakes  the  first 
the  velocity  immediately  to  the  right  of  the  wave  front  is  v.  The 
velocity  of  the  entire  string  to  the  right  of  the  wave  cannot  change 
abruptly  from  V2  to  v;  hence,  unless  these  quantities  are  equal, 
the  merging  of  the  two  waves  does  not  take  place  smoothly,  and  a 
tightening  or  slackening  of  the  string  occurs.  The  corresponding 
change  in  stress,  and  hence  strain,  constitutes  a  disturbance  which 
causes  a  corresponding  disturbance  to  be  propagated  along  the 
string  to  the  right  away  from  the  wave  front. 


Substituting  (47)  and  (49)  in  (46),  and  (51)  in 

(  52) ,  we  obtain 

AT  - 

nrz  z 

Since  the  stress-strain  curve  is  concave  upward 

it  is  evident  that 

, 

or 

(Ta.  €  ,  >  <T  j  €  2  . 

It  follows  that 

>  \FT7j. 

0  >  -  +  a.  V<r, 

<s «•«;«, -rte,-*i iC2 

-V5S)‘  > 
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Fig,  12*  Stress- strain  curve  showing  chords  that  determine  velocities  of  vertical 
wave  fronts* 


o  Lx 

Fig.  13.  Wave  front  of  arbitrary  but  constant  shape. 
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In  view  of  (53)  this  may  be  written 

flT  7*  •  (55) 

It  follows,  as  described  on  page  33,  that  the  merging  of  the  two 
waves  is  not  smooth,  but  causes  a  disturbance  which  travels 
along  the  string  in  the  opposite  direction  from  that  of  the  wave. 

Note,  The  situation  portrayed  in  Fig.  11  can  be  achieved  by 
using  a  stress-strain  curve  composed  of  two  straight  lines. 

One  passes  through  the  origin,  and  is  valid  for  0  <  e  <  e1;  the 
other  is  steeper,  and  is  valid  for  E\  <  £  <  £z. 

Effect  of  a  Change  of  Shape  of  the  Strain- Wave  Front,  From  Fig.  12 
and  the  equations  (47)  ,  ( 49) ,  (  51) ,  and  (55)  we  saw  that  when  one 
abrupt  strain  wave  front  catches  up  with  another  there  is  a  sudden 
slackening  of  the  string  behind  the  wave  front,  and  a  sudden  increase 
in  the  velocity,  or  surging  forward  of  the  wave.  This  suggests 
that  the  wave  front  loses  its  abruptness  and  constancy  of  velocity. 
Accordingly  let  us  consider  the  case  of  a  wave  front  of  constant 
shape  which  moves  in  translation  with  a  velocity  c(t),  no  longer 
required  to  be  constant,  Fig.  13. 

At  time  t  =  0  let  the  displacement  s  (x,  o)  be 
denoted  by  f ( x) ,  thus 

rfC*).  (56) 

At  present  we  shall  leave  f(x)  arbitrary.  At  time  t  this  disifa.ce- 
ment-position  curve  moves  a  distance  C  c  dt  to  the  left;  hence 


(57) 


In  Appendix  2  an  alternative  point  of  view  is  presented  (but  not 
worked  out)  in  which  the  sudden  increase  in  propagation  velocity 
does  not  occur. 
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(58) 


v  =j|  (*■■*■/„*•&)  =,c€- 


At  time  t  let  the  front  extend  from 
x  =  0  to  x  =  L;  then  the  momentum  of  that  part  of  the  string 
which  composes  the  front  is 

M=f  (6o) 

where  A  is  the  cross-section  area  of  the  string,  and  s^  is  the 
displacement  at  the  right  hand  boundary  of  the  front.  The 
stress  at  this  boundary  is 


<r l.  - 


(61) 


where  £ 


L 


is  the  strain  at  this  boundary. 


Applying  the  momentum  theorem  (37)  to  that 
region  which  contains  the  front,  and  which  moves  with  the  front, 
we  ohtain 


^(Af  caJ  +A[x(c  $l)  -  A 


or,  canceling  A  and  noting  that  all  qiantities.  except  c  and  t  are 
constant. 


(62) 
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This  differential  equation  can  be  solved  exactly;  however,  it 
is  evident  by  inspection  that,regardless  of  its  initial  value,  c 
is  going  to  approach  the  value  obtained  by  placing  dc/dt  m  0 
in  this  equation.  This  final  value  of  c  is 


regardless  of  the  shape  of  the  front. 


(63) 


We  note  that  this  value  agrees  with  (47),  which  applies  to  an 
abrupt  wave  front. 


Next,  let  us  suppose  that  c  has  attained  its 
final  value  ( 63).  We  shall  consider  the  energy  balance  of  the 
same  region  to  which  we  applied  the  momentum  theorem,  namely 
that  which  contains  the  front  and  moves  with  the  front.  Since  the 
region  moves  with  the  front,  we  see  that  the  rate  of  increase  of 
kinetic  energy  and  the  rate  of  increase  of  potential  (strain)  energy 
in  this  region  are  both  zero.  Also,  noting  (59)  we  see  that 

External  power  input  =  A  G~]_  ^  (  64) 


Kinetic  energy  loss  per  second 
through  convection  = 


n<rLceL 


due  to  (63).  Also, 


I*^L. 

-/C  Aj  <rte 


(65) 


~  A  G~cuot  & 


(66) 
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For  there  to  be  an  energy  balance  it  is  necessary  and 
sufficient  that 


=  -s'  r. 


(67) 


Here  <r  is  the  average  value  of  <r  ,  averaged  on  c 
from  0  to  eeL»  We  note  that  the  condition  (67)  is 
independent  of  the  shape  of  the  front.  In  fact  the  above 
relations  agree  with  those  on  pages  29  and  30,  which 
pertain  to  an  abrupt  wave  front.  It  follows  that  the 
situation  in  regard  to  energy  balance  cannot  be  improved 
or  even  altered  by  changing  the  shape  of  the  front. 

Overtaking  of  a  Wave  Front  by  an  Infinitesimal  Wave,  On 
pages  31  to  35  it  was  shown  that  with  a  concave  upward 
stress- strain  curve  one  abrupt  front  overtaking  another 
will  not  smoothly  merge  with  it.  To  be  certain  that  the 
situation  doesn't  appear  different  when  the  overtaking  front 
is  treated  as  a  series  of  infinitesmals  let  us  consider  the 
case  where  an  infinitesimal  wave  comprising  a  strain 
increase  de  and  moving  with  a  velocity 


/C,  >  /C  (68) 

overtakes  a  wave  front  comprising  a  strain  increase  € 
and  moving  with  a  velocity  c,  as  indicated  in  Fig,  14, 

From  (38),  (41),  and  equation  (138)  of  the  appendix  we 
have 


or  -  / c-  £  ; 


(69) 


JL.  -  >  (70) 

,C,=  (71) 
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Fig.  14.  Vertical  wave  front  being  overtaken  by  an  infinitesimal  wave. 
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ft 

where  <r  =  <r  (e),  and<r(£)>0.  The  stress-strain  curve  is 
thus  concave  upward,  and  c  increases  with  £.  Differentiat¬ 
ing  (69),  <r  ,  and  (70),  we  obtain 


—  /C  cL£  +■  £  etc* j  (72) 


JL<r  -  <r' (0  oie  , 


dc  - 


j.  ct<r 

«*■ 


or,  noting  (70)  and  (71), 


(73) 


(74) 


Equations  (72)  and  (74)  give  the  increases  in  v  and  c, 
respectively,  when  the  infinitesimal  wave  catches  up  with 
the  main  wave,  thereby  increasing  it  by  d£.  Substituting 
(74)  in  (72)  now  gives 

JUr--  4  c'~£-  di 

2. 


i/C 


(75) 


Since  Cj  >c,  which  follows  from  (70  ) , 
(71) ,  and  the  fact  that  the  stress-strain  cvirve  is  concave 
upward,  we  see  from  (74)  and  (75)  that  dc  and  dv  are  both 
positive.  These  results  are  in  agreement  with  those 
obtained  on  page  s  3  2  and  35. 
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Consideration  of  a  Discrete  System  of  Masses  and  Springs 
as  an  Approximation  to  a  Continuous  System.  We  have  shown 
(see  pages35r38  )  that  no  wave  front  of  any  constant  shape 
can  exist  when  the  stress-strain  curve  is  concave  upward. 

In  order  to  get  some  idea  of  what  is  going  on,  without  becoming 
involved  with  a  lot  of  calculations,  we  shall  consider  the 
discrete  system  of  masses  and  springs  shown  in  Fig.  15.  The 
masses  are  represented  by  the  numbered  circles,  are  all  equal, 
and  extend  indefinitely  to  the  left.  The  nonlinear  springs  which 
connect  adjacent  masses  are  all  alike,  and  are  represented  as 
sagging  cords.  Until  the  slack  is  taken  up  these  cords  carry  no 
tension  whatsoever;  but  as  soon  as  the  slack  is  taken  up  they 
behave  as  infinitely  stiff  elastic  springs.  The  corresponding 
stress-strain  curve  is  shown  in  Fig.  16,  though  here  the  role 
of  c  is  played  by  tension,  and  the  role  of  e  is  played  by  the  loss 
of  slack,  or  relative  displacement.  The  masses  are  all  initially 
at  rest  as  shown;  and  starting  at  time  t  =  0  the  right  end  moves 
to  the  right  with  a  constant  velocity  vQ.  We  note  that  the  stress- 
strain  curve,  though  consisting  of  a  broken  line,  lies  in  the 
category  of  concave  upward. 

Whenever  the  slack  is  taken  up  in  any  of  the 
cords,  impact  conditions  prevail.  If  the  cord  is  that  to  the  right 
of  mass  1,  this  mass  changes  its  velocity  abruptly  from  0  to 
2  v0.  If  the  cord  lies  between  two  masses,  those  masses  abruptly 
exchange  their  velocities.  These  facts  are  enough  to  enable  one 
to  draw  the  displacement- time  curves  for  the  various  masses,  as 
shown  in  Fig.  17.  In  drawing  these  curves  we  work  from  left  to 
right,  and  note  that  slope  is  velocity.  The  slack  is  four  divisions; 
and  when  the  vertical  distance  between  adjacent  curves  attains 
this  value,  impact  occurs  and  slopes  change  abruptly. 

As  time  goes  on,  the  number  of  masses  involved 

in  the  motion  increases;  also,  since  slope  is  velocity  we  see  that 

of  these  masses  alternate  ones  are  at  rest,  and  the  other  alternate 

ones  are  moving  with  a  velocity  2vq,  this  being  true  at  any  time. 

The  velocity  distribution  along  the  masses  may  hence  be  regarded 

as  composed  of  a  constant  velocity  vq  plus  one  which  alternates 

between  v  and  -v  as  we  pass  from  mass  to  mass.  The  alternat- 
o  o 

ing  velocity  does  not  contribute  to  the  overall  motion,  but  it  doubles 
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Fig.  17.  Displacement-time  curve  for  string  of  Fig.  15  subjected  to 
constant-velocity  impact. 


the  kinetic  energy  contained  in  those  masses  which  are  involved 
in  the  motion.  We  see  that  although  the  disturbance  moves  to 
the  left  at  a  uniform  rate,  it  leaves  behind  it  not  a  series  of 
masses  all  moving  with  the  same  velocity  vQ,  but  a  motion  having 
in  addition  a  violent  oscillation  of  amplitude  vQ.  We  have  thus 
encountered  a  situation  similar  to  that  which  we  have  been  led  to 
believe  takes  place  in  a  string  whose  stress-strain  curve  is  con¬ 
cave  upward,  after  a  shock  wave  has  formed. 

Figure  18  is  similar  to  Fig.  17,  but  corresponds 
to  a  case  wherein  the  velocity  of  the  moving  end  is  altered  during 
the  course  of  the  motion.  It  was  thought  that  possibly  such  an 
alteration  would  cause  a  permanent  scrambling  of  the  pattern  of 
motion;  but  such  is  evidently  not  the  case. 

In  an  attempt  to  approach  the  continuous  distribu¬ 
tion  let  us  now  divide  each  mass  into  n  equal  masses,  and  dis¬ 
tribute  these  evenly  along  the  string,  the  total  mass  per  unit  length 
remaining  the  same.  The  slack  in  the  string  between  adjacent 
masses  is  thereby  also  reduced  in  the  ratio  l/n.  Proceeding  as 
before,  we  again  obtain  a  motion  wherein  alternate  masses  of  those 
involved  in  the  disturbance  are  at  rest,  and  the  remaining  alternate 
masses  are  moving  with  a  velocity  2vQ.  Impacts  occur  n  times  as 
frequently  as  before,  and  the  disturbance  progresses  to  the  left  at 
the  same  constant  rate  as  before.  In  the  disturbed  region  the 
velocity  distribution  may  be  considered  to  be  composed  of  the 
constant  velocity  vQ  plus  a  violent  oscillation  of  amplitude  vQ  as 
before;  however,  the  wave  length  is  now  l/n  times  what  it  was 
before.  If,  now,  n  is  increased  indefinitely  in  order  to  approach 
a  continuous  mass  distribution,  this  wave  length  approaches  zero; 
and  the  oscillation  retains  its  violence,  but  becomes  infinitely  fine. 

It  may  be  mentioned  that  in  the  case  of  discrete 
masses  even  if  the  springs  were  linear  an  end  disturbance  would 
not  be  propagated  without  distortion;  for  we  are  dealing  with  the 
mechanical  analogue  of  a  low  pass  filter.  Different  frequencies 
are  propagated  with  different  velocities;  and  above  a  certain 
cutoff  frequency  are  not  propagated  at  all. 
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More  explicitly  let  us  consider  the  linear 
system  of  Fig.  19.  The  masses  are  all  of  magnitude  M,  and 
the  spring  constants  are  all  of  value  K.  The  masses  are 
numbered  from  left  to  right,  and  the  displacements  toward  the 
right  are  indicated  by  sj.  The  equations  of  motion  are 

.  <76 

For  vibratory  motion  we  place 


A:  - 


(6 o~t  +Y) 


(77) 


hence,  (76)  becomes 

where 


(78) 


(79) 


We  note  that  u>  is  the  resonant  angular  frequency  of  any 
mass  if  the  two  adjacent  masses  are  held  fixed.  If 


<  I 


the  solution  of  the  difference  equation  (78)  is 


(80) 


where 


t 

-  C  i  sCAXL  -Cot.  -f 


Atsvu  cK 


-I 

-  Aurtb 


0 


(81) 


(82) 
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Equation  (81)  in  (77)  gives  a  standing  wave  of  arbitrary- 
amplitude,  phase,  and  position;  however,  by  superimposing 
two  of  these,  thus 


~  -C  CsHCot  -  A.Co4>  X.C co'h 


-  ,C  (  J.  c*.  -  oo£)  > 


(83) 


we  obtain  a  traveling  wave.  This  wave  moves  to  the  right 
along  the  i-axis  with  a  velocity 


OJ 

^  *  “3T  > 


(84) 


which,  in  view  of  (82),  depends  upon  the  angular  frequency  0). 
In  this  wave  only  integral  values  of  i  pertain  to  the  masses, 
and  thus  have  direct  physical  significance.  Since  co  is 
positive,  ( 80)  is  equivalent  to 

0  <  OJ  <  yfi  .  (85) 

If 

to  y  ifz  lo  c 

( 80)  is  violated;  and  in  place  of  the  traveling  wave  (83)  we 
obtain  the  attenuated  wave 

•  * 

/C('0  £  (a)t  4  y)  (87) 

where 

o(  - (.60^  0* 
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Finally,  if  CO  =  0  there  is  no  motion.  We 
now  see  that  a  wave  is  propagated  only  if  (0  lies  below  the 
cutoff  value  C0q,  in  which  case  the  velocity  of  propagation 
depends  on  CO. 

Oscillations  in  the  Continuous  System.  The  oscillations  which 
we  have  encountered  in  the  discrete  system  just  considered 
(see  page  45)  suggest  a  possible  way  out  of  the  difficulty 
which  we  found  in  connection  with  the  lack  of  energy  balance 
in  the  propagation  of  an  abrupt,  or  square  wave  (see  page 31). 

It  may  be  that  the  extra  energy  which  is  available  in  the  case 
where  the  stress-strain  curve  is  concave  upward  is  used  to  pro¬ 
duce  oscillations  in  that  part  of  the  string  which  has  been  passed 
over  by  the  wave  front.  In  investigating  this  possibility  we  shall 
consider  the  string  in  the  unstretched  state  as  extending  from 
x  =  -  00  to  x  =  0;  and  at  time  t  =  0  the  right  end  of  the  string  is 
suddenly  given  a  constant  velocity  vQ  to  the  right. 

We  shall  consider  a  resulting  strain  wave  of 

the  form 

£  ?  ^  ^  "f  (2*j  AXMj  H  (7C  (X 

(89) 

where  eQ»  cq,  _Q_  ,  and  the  a*s  are  constants,  and  where 

iir 

_Q.  -  -y-  ,  (90) 

\  being  the  wave  length.  In  addition  to  the  constant  value 
e  we  are  thus  providing  for  the  possibility  of  two  sinusoidal 
waves  of  wave  length X.  moving  in  opposite  directions  with 
velocity  cq.  If  these  two  waves  were  infinitesimal,  we  know 
that  cQ  would  be  the  velocity  corresponding  to  a  strain 
(see  appendix  1,  page  62).  In  any  case  we  should  expect  cq 
to  approximate  this  velocity,  which  is  -/cr(  £0)/ f.  We  shall 
try  to  determine  the  various  constants  in  (89)  so  as  to  satisfy 
various  dynamical  conditions,  one  of  which  is  that  the  energy 
unbalance  which  we  encountered  on  page  31  be  removed. 
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The  differential  equation  which  determines 


the  motion  is 


r 


^  <r 
s  % 


or 


p  =  a-Vo  — 
f  Tit1  1  J  - 


At  x  =  0  we  have 

Aj  -  ar0£  a^eL  -  O  x  -  O^* 

hence  (92)  gives 


a  e 


—  £)  9^  — •  0  • 


(91) 


(92) 


(93) 

(94) 


Substituting  (89)  in  (94),  we  obtain 

JX  Co-i 

'*a'3  CoiJliCiiL  ~~  4*+  SL  C0~k  J  —  0  . 

This  relation  will  be  satisfied  for  all  values  of  t  if  and  only 
if 

£L3  -  ~  ^ I  )  CL±-  <X>2_  > 


(95) 


due  to  which  (89)  becomes 

£  -  (l0  +  a([^c-flC^"c-wi)-^-n-C)t+A«^)] 


4- 
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This  may  be  written 

€  =  -  Xa.tCMijSi7.4i^n.c0t  +&<*.*  co4Sit..± 

=  €«,  +XC0<uSl%  (-<LtM^SlCot  +CL-tc^sx.ctib),  (9 7) 

or 

£  -  £o  -f-  -TL^t  Am**-*  (/1c6£  +■  ^ )  (98) 


where  and  y  are  constants.  We  see  that  e  is  composed 
of  a  constant  and  a  standing  wave  of  wave  length  \  = 
and  angular  frequency  0)  =Jlc0. 


Integrating  (98)  with  respect  to  x  holding 
t  constant,  we  obtain 

A,  -  £0%  +  jt*^SL7:AA^j(£).c0±  +  y)  i-'+W  (99) 

where  i]i  (  t )  is  an  unknown  function  of  t#  But 


s  =  v  t  when  x  =  0. 
o 


(100) 


which  with  ( 99)  gives 


4*  (t)  =  vQt; 


hence 


-  C  ~Y  -h  — “  -4 Xvu  _fl”X  +y)  i 

'  co  A  Jt- 


Differentiating  partially  with  respect  to  t,  we  obtain  the 
following  expression  for  the  velocity 

ax  =  rif'  +  6,/c0^  At:  c** 

The  kinetic  energy  per  unit  length  of 


string  is 


(101) 


(102) 


■A  r  -X |^ATo  +  i/!T#  €- j  £,c  coo  (-^-^0^ 

•f  ,t<nx(£lCo±  +5f)  J  (103) 


where  A  is  the  cross-section  area. 
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The  average  of  this  expression  taken  over 
\  =  2tr/fi  is  y  +  xjt 


a  wave  length 


=  ££[ail  +  ^rc«.1  (■**.*  +*)]; 


(104) 


the  average  taken  over  a  time  period  2tt/  cQ  is 

six,,  ' * '*£• 

p  '  att, 

£t> 


-I 


4- 


^ l&a  .  ■  2 

si  7C 


<■]> 


(105) 


and  the  average  over  both  is 

-r  jjjSi*  /*•* 

-  +5FV 

,  IA.  [/v;2  +  -A_|f  j .  (106) 

Here  and  tQ  are  arbitrary  values  of  x  and  t,  respectively. 

The  strain  energy,  or  potential  energy,  per 
unit  length  of  string  is 


iT 


T  £fi- 

J  ~ 

£  a 


6 

<r(i 6)  dt 


(107) 


hence,  before  this  can  be  averaged  over  a  space  period  X.  or 
a  time  period  2TT /fLc^cr  (  e)  must  be  specified.  For  the 
present  let  us  place 


cr(e)  -  £0£  +  E,€z 
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(108) 


where  E  and  Ej  are  both  positive  constants.  We  note  that 
this  stress-strain  curve  is  concave  upward.  Equation  (108) 
in  ( 107)  now  gives 


(109) 


Substituting  (98)  in  this  expression,  we  obtain 

~ ~°( +y) 

+  6  2  'Ctrtf'JITC  +#)] 

■f  — ~ +3^0  (Ac*t 

+ 1  C-A-C'i  +  y)  (no) 


The  average  of  this  expression  over  a  wave  length  is 


h--£rf 


J2-  I'*** 

*o 


§U 

~  [V  +  u(-n.Cob  +  y)] 

=  -fcA6^Eo-^E/#0) 

+  *AS(\£0+2£,60)^l(-n.£o;fc  (HI) 


the  average  over  the  time  period  a/IT/riAiis 

r%  +■ 

ir&I  i* 


^  -^A £,  [ft,  /c^ru 

+i  A  +1  £<  ^ 


(112) 
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and  the  average  over  both  is 

*  ,  liT  ^  *42. 

-A-  /*£<>  * 

^>pr~T7 r-j^  ^  §_<tt  d\ 


—  (^o 

+  7  A  6,  (5C  +i.E,  ^0)  ,  (113) 


Adding  (104)  and  (111)  we  obtain  the  following 
expression  for  the  average  of  the  total  energy  taken  over  a 
wave  length. 

ta*  +i**l(3£°f2£>e°) 

n  *€»*(£•*  *£i  <o  *  Li  r*  CcJ 

-1^6,  (Ee+AE/djJ/xa  2(JXC0±  4  y).(114) 


There  is  no  mechanism  for  transferring  energy  to  and  from 
an  appreciable  length  of  the  string  at  a  frequency  which 

is  twice  the  frequency  of  vibration  of  the  standing  wave.  In 
order  to  avoid  a  situation  where  this  is  implied,  we  equate  the 
coefficient  of  the  cosine  in  ( 114)  to  zero,  thus 


=  O  <115> 


or 


£*  +2£,ec 


(116) 
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This  expression  determines  c  .  Noting  (  108)  we  see  that 
(116)  may  be  written 

>c0  =  y<T(fo 7} ,  (in) 

which  is  the  wave  velocity  corresponding  to  a  strain  eQ. 

This  is  in  line  with  the  conjecture  made  on  page  49.  In  view 
of  (115),  (114)  reduces  to 

=if  + 

+  ifi£(3Bo  *■*£,£.)+ (ns) 

or,  substituting  (116)  for  cq, 

T*+^»=i  f  A  j 

A  6|  (E-o  j 

Adding  ( 105)  and  (112)  we  obtain  the 
following  expression  for  the  average  of  the  total  energy 
per  unit  length  taken  over  a  time  period  2  IT  / SL  /Ce  . 

V^P  "ifA 

(120) 

—fcA£ii'(E0+Z£,6o')]  * 

But  here  the  coefficient  of  the  cosine  is  zero  due  to  ( 115); 
hence,  (120)  duplicates  (118),  or  (119).  We  thus  have 
finally,  noting  (106)  and  (113) 

— if 4  ”5  ^  ^  ®  (3  £o^  i2i) 
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Next,  let  us  suppose  that  the  disturbance 
moves  to  the  left  with  an  average  velocity  c.  We  do  not 
require  that  there  be  an  abrupt  wave  front,  nor  that  the  wave 
velocity  c  be  constant.  In  fact  at  this  point  we  do  not  have  to 
know  the  detailed  nature  of  this  front  -  only  that  there  is  a 
transition  region  of  finite  length  between  the  unstressed  string 
and  the  string  in  the  stressed  state  described  above  ( see  eq  98), 
which  region  moves  to  the  left  with  an  average  velocity  C.  We 
shall  now  proceed  to  obtain  three  conditions  which  must  be 
satisfied  for  the  conservation  of  mass,  momentum,  and  energy. 

Conservation  of  Mass.  The  total  stretch  which  occurs  in  the 
transition  region,  or  "  front"  ,  will  be  denoted  by  s^,.  We  do 
not  know  this  quantity;  however,  we  do  know  that  it  is  bounded. 
Referring  to  Fig,  20,  let  x^,  be  the  abscissa  of  the  right  hand 
boundary  of  this  front;  then  ( 101)  gives 


.  (122) 

But 

~  +  (123) 

where  xx(t)  is  bounded.  This  in  (122)  gives 

-fl-7Cp  (-fl  C.0  jt,  -f  y)  t  (124) 

As  time  goes  on,  all  terms  in  this  equation  are  bounded 
except  the  first  and  last  terms  on  the  right  hand  side. 

These  terms  must  cancel  each  other;  otherwise  the 
equation  will  ultimately  be  violated.  It  follows  that 
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(125) 


C  =  —  • 

Conservation  of  Momentum.  The  total  momentum  in  the 
string  is 


M  --  M, 


'r 


AT  dsf-  , 


(126) 


where  is  the  momentum  in  the  front,  and  is  known  only 
to  be  bounded.  Substituting  ( 102)  in  this  expression  we  obtain 

o 

M=MF  +ApJ[nre>i6,fi0AU>SLX 

=.  Mp  +  Af  1 

+  ( 12?^ 

This,  however,  must  be  equal  to  the  total  impulse  applied 
at  the  right  hand  end;  hence 


M-- a/ 


(128) 


Placing  x  =  0  in  (98),  then  substituting  in  ( 108)  and  ( 128), 
and  then  adding  an  unknown  but  bounded  term  to  allow 

for  the  fact  that  a  short  time  may  be  required  for  £  to 
attain  the  form  (98),  we  obtain 

A. 

Ji  =  A l  {£*[*.  +  +*)]  +  E,[tZ 

2(a£.,£ ( 1 29 ) 
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Replacing  x^,  in  ( 127)  by  ( 123),  then  equating  the  unbounded 
terms  in  (127)  and  (129),  we  obtain 

=  A  L£.  +-r)j 

-  E-C^e  (130) 


Conservation  of  Energy.  The  energy  imparted  to  the  string 
by  the  force  at  the  right  end  is 


± 


(131) 


This  expression  differs  from  ( 128)  only  by  a  factor  of  vq, 
and  can  hence  be  obtained  by  multiplying  (129)  by  vQ, 
thus 

£  -  A nrQ  {fcK  4  E, Cf*  +  )& 

~/C*a  *] 

The  total  energy,  kinetic  and  potential, 
contained  in  the  string  is 


T  +  §^(V?A)2-*  +  £b 


(133) 


where  <5 „  is  bounded.  Substituting  (119)  in  (133),  then 
equating  the  unbounded  terms  of  (132)  and  (133),  we 


obtain 

+iAc‘(3E.+J.E,0 

+  ±A«,*(E.AA-eiE»)]£  =AW„[_E„fc<. 

*■  J  cx  1 

+•  £/(^c  +  "1  )J) 
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(134) 


or  noting  ( 130) 

■Z$(3E.*-12lQ+i<Xe.+lSl6.)  =  i  frfiS 
i£‘(3^i£,Q  =  i  f  V,' . 

Substituting  (125)  in  (130),  we  have 

A  +  F,  f .)  +  6/ *f ,  €„  =  1  f  ^ , 

or 

€,1  = 

Substituting  (136)  in  (134),  we  obtain 

finally 

E,6?(3£***£i6o)  +3(£<>  *af,60) 

•  L  (£0-t  £i£C>')]  =  ?  f1  , 

or 

4  £*  £+ (*£,£»** 

~  3  PWoE,  eo  -3(°^0X  £e  -  C> . 

1  1  (137) 

This  quartic  equation  determines  eQ,  after  which  £j  is 
given  by  ( 136) ,  and  5  by  (125). 

Thus  we  have  shown  that  for  a  particular 
concave-upward  stress-strain  curve  as  given  by  (108)  it 
is  possible  to  specify  the  constants  in  an  oscillatory 
solution  (98)  to  a  tensile  impact  problem  so  as  to  give  an 
energy  balance.  The  implication  is  that  for  any  concave- 
upward  stress- strain  curve  the  solution  would  be 


(135) 


(136) 
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oscillatory.  However,  since  we  have  seen  that  the 
iterative  procedure  used  in  our  computer  calculations 
will  produce  oscillations  at  a  discontinuity  even  when  they 
are  not  required  for  energy  balance  it  would  be  pointless 
to  try  to  account  for  the  oscillations  observed  in  the  com¬ 
puter  results  from  energy  considerations  alone. 

Before  leaving  this  theoretical  analysis 
of  oscillations  it  shoiil  be  pointed  out  that  in  some  instances 
problems  have  been  left  without  working  out  all  the  details 
and  without  utilizing  all  the  known  conditions  of  the  problem. 
This  was  done  because  we  were  anxious  to  get  on  with  solutions 
involving  creep  of  the  material.  We  believe  that  creep  is 
always  present  and  that  a  time-independent  model  can  never 
accurately  represent  reality.  Hence,  although  a  more  complete 
study  of  the  time- independent  model  could  be  made,  and 
although  such  a  study  might  help  in  the  understanding  of  the 
propagation  of  actual  strain  waves,  we  will  not  develop  the 
time- independent  analysis  further  at  this  time. 

8.  Conclusion 


It  has  been  shown  in  the  present  report  that 
a  non-linear  time- independent  tensile  impact  problem  may  be 
solved  by  an  iterative  computer  procedure  and  that  the  result 
agrees  with  the  solution  obtained  by  the  method-of- characteristics. 
In  either  method  difficulty  may  arise  if  the  stress-strain  curve 
becomes  concave  upward.  In  the  method-of-characteristics  these 
difficulties  are  overcome  by  assuming  shock  conditions  (and  by 
implication  at  least  assuming  the  production  of  heat).  In  the 
iterative  computer  procedure  the  difficulties  cause  oscillations. 
Although  partly  computational  it  has  been  shown  that  in  many 
ways  these  oscillations  are  a  more  correct  solution  than  the 
assumed  shock. 


We  are  now  ready  to  apply  the  iterative 
computer  procedure  to  the  time  dependent  non-linear  case  which 
we  believe  is  closer  to  reality  than  anything  discussed  so  far. 
This  will  be  done  in  Report  m. 
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Independence  of  Infinitesimal  Strain  Waves,  The  purpose  of 
of  this  section  is  to  check  the  validity  of  the  idea  that  in  the 
case  where  cr  =  <r(  e)  each  infinitesimal  strain  increment 
travels  with  its  own  velocity 

c(e)  =  y<r'(0/f 

independent  of  the  others.  This  idea  is  suggested  by  the 
differential  equation 

_  a  <T 

P  z t1  ex 

written  in  the  form 


~  l 

z 

Tx1  (141) 

where  c  is  given  by  (138).  For  an  infinitesimal 
disturbance,  c  is  essentially  constant,  (141)  becomes 
the  wave  equation  closely,  and  it  appears  that  the 
disturbance  is  propagated  undistorted  at  velocity  c  (c). 

If  the  disturbance  is  an  infinitesimal  strain  increment  de, 
we  have  the  above-mentioned  idea.  Although  the  above 
provides  a  suggestion,  it  leaves  a  great  deal  to  be  desired 
by  way  of  proof.  In  what  follows  we  shall  suppose  that  the 
idea  is  valid,  and  shall  determine  whether  or  not  the  strain 
wave  which  it  gives  satisfies  the  differential  equation. 


(140) 


(139) 


(138) 


At  time  t  let  the  strain  wave  be  represented 
by  the  curve  in  Fig.  21.  As  time  goes  on,  the  future  shape 
of  this  curve  is  determined  by  the  motion  of  the  individual 
strain  increments  —  each  shearing  by  the  others  at  its  own 
velocity  c(£),  as  indicated.  In  time  dt  each  strain  incre¬ 
ment  de  moves  to  the  left  a  distance  c  dt.  It  follows  that  at 
abscissa  x  (of  the  unstretched  string)  the  strain  increases 
during  this  time  by  the  amount 

cU  ^  ||  *-(0  oik  ; 

hence,  dividing  by  dt, 

<142> 

which  is  a  quasilinear  first  order  partial  differential  equation. 
Equation  (142)  does  not  resemble  the  desired  equation  (  140) 
or  ( 141) ;  and  cannot  be  transformed  into  it  in  any  evident 
manner.  However,  there  is  still  the  possibility  that  the 
solution  of  (142)  will  satisfy  (141).  We  shall  therefore  solve 
(  142).  Equation  (  142)  may  be  written  in  the  standard  form 

-/C(£)  -  O  ,  (143) 

from  which  we  see  that  the  characteristics  are  given  by  the 
equations 

Jjt  (144) 

T~  =  -*(0  ’  0 

These  give 

£  -r  /C-2.  -  }  (145) 
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Fig.  21.  Strain-wave  front  composed  of  incremental  strains  each  moving  at 
its  own  velocity. 
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whence 

^TC-f/O  (£*.)  ebb  -  o 

7(  -+C(Ct)£  -  C,  -  zmu^LuJ: .  (146) 

Replacing  C2  by  e  we  obtain  for  the  general 
solution  of  ( 143),  and  hence  ( 142) 

£  -  -f  [X  4  <£(£)£]  (147) 

where  f  is  an  arbitrary  function.  Initially  t  =  o  and 
(147)  becomes 


6  =  fOO  (148) 

from  which  we  see  that  f(x)  is  the  arbitrarily  specified 
initial  distribution  of  e. 

From  (147)  we  obtain 

JL6  =  (fOCrf*  +  /C* *  ^  +  ^ett)  (149) 

where 

r  '  -  {' 1%  + 

^  =  £.£€) 

Equation  (149)  may  be  written 

/ 


from  which  it  follows  that 
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a  6  _ 

f 

(152) 

*7 c 

l-f'A*  ’ 

3  6  _ 

,cf' 

(153) 

l-fWst  * 

We  note  that  these  quantities  satisfy  the  differential 
equation  (142).  In  order  to  check  the  differential  equation 
( 141)  we  need  S^s/l  t^hich  we  shall  now  proceed  to  calculate. 


Since 

,  - 

£  ~  37.  (154) 

we  obtain  s  by  integrating  e  with  respect  to  x  holding  t  constant, 
thus 


-A  =  f f  [l(+ JjC  +-  ^4j(0,  t)  } 


(155) 


where  s(o,t)  is  the  value  of  s  at  x  =  o  and  time  t.  It  follows 
that 


f -4  -t'C0)*)' 


Since,  noting  (152)  and  (153) 

a 

ax 


( 156)  may  be  written 

aA 
air 


rx 

~Jf  [~X-  +'£'($)  t-M)  5x^*  ^ 


+aj-4'(6'jfc). 


(156) 


(157) 


(158) 
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Integrating  by  parts,  we  obtain 


r 

-  (-f  Ctc  +  UQtl  *C0j[ 


Differentiating  again,  we  have 

=  j-f  [X  4C&)  £]  «&>$£ 

-f  f  [7C  +fi (*)&]<.(£)[£■(£) 

-f-  J^*z.  ^  (  0  )  ^0 

'[%  +  <(<>)£ ]  fcCO  %|  ]^)|| 

4-f  W+cCOi  J  4fo)S  ||  4f£x*c(f)^)|^jcL*.j 

Here  in  the  first  bracket  we  note  from  ( 147)  that 


f!  f -«*]  im +*a'«)  H  j  ; 
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also,  since,  noting  (152)  and  (153), 

ii 


it  follows  that 


jU. 


C  ^  ^  (0X)‘ 


(159) 


(160) 


(161) 


(162) 


(163) 
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Furthermore 

£  (€)  37C  j 

2  f  tx  +-C(€)1]  f  I  +  *AVO^ 

4-  f£X  +C&)*]  ^  ^  ^Xx  * 


In  view  of  these  relations  the  integrand  of  the  integral  i 
( 1 60)  becomes 


i-fCX +<'(*)£]{  fiSC*)  14] 


and  ( 160)  becomes  ^ 

=  5  -ft**  'C&t] &'(*)§%■ 

*£.^c°,±) 

_  ^ <c«)  5x/ f  (*+£•(()£]<'  (*)  _/ 

*f£**'C/0*j/*,VO  7 

Integrating  by  parts,  we  now  obtain 

4^4  ^+[x-k(6)£] §i 

x-tcc^tjc'co 

+  £m°,±) 

+~J  j’etojj  f[X4x60fc] 

-  f  (*+*<«)*]  [■c-'CfVyf]  ix .  I 


(164) 


(165) 


(166) 


(167) 
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Noting  ( 162)  we  see  that  the  first  and  last  terms  in  the  first 
bracket  cancel  each  other;  also,  the  two  terms  in  the  integrand 
cancel.  We  then  have 


— 

" 


£ 


Mo ,i)) 


(168) 


which  in  view  of  (162)  becomes 


(169) 


or 


pc=o 


(170) 


We  saw  from  (147)  that  the  process 
indicated  in  Fig.  21  is  sufficient  to  determine  e  if  f  (x), 
the  initial  distribution  of  e,  is  given.  From  (155),  however, 
it  is  evident  that  s  is  determined  only  to  within  an  arbitrary 
function  s  (o,t).  If  for  x  =  o  we  choose  for  this  function  one 
that  is  physically  possible,  then  the  differential  equation  (141) 
is  satisfied  at  x  =  o,  the  unknown  left  side  )  _ 

being  determined  by  the  known  right  side 

o  ; .  • 
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It  follows  that  the  large  bracket  in  (170)  vanishes,  and  ( 170) 
becomes 


<D  X* 


) 


which  duplicates  the  desired  differential  equation  (141),  We 
have  thus  shown  that  the  idea  of  independent  propagation  of 
infinitesimal  strain  increments,  as  described  on  page  62  and 
indicated  in  Fig,  21  is  valid. 

In  Fig,  21  f  (x)  is  shown  with  a  positive  slope. 
This  slope  may,  however,  be  negative,  as  shown  in  Fig.  2 £— 
the  above  mathematical  treatment  is  valid  in  both  cases. 

From  the  denominators  of  (152)  and  (153)  we 
see  that  if  the  slope  of  f(x)  is  positive,  and  if  c*  (c  )  is 
positive;  then  after  a  certain  time  has  elapsed  defix  and 
ae/  at  become  infinite.  For  strain  £  this  time  is  given  by 

/  -  o 

or 


£  - 


(171) 


Here  f*  is  the  initial  value  of  Be/  <kfor  the  specified  value  of  £, 
as  is  evident  from  ( 150).  The  smallest  value  of  t  given  by 
(171)  for  our  range  of  values  of  e  is  the  time  required  for  a 
singularity  to  develop.  At  this  time  the  validity  of  the  above 
mathematical  treatment  ceases;  and  at  subsequent  times  some 
extension  of  this  treatment  is  necessary. 
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Appendix  2 

Here  some  additional  comments  will  be 
made  on  the  problems  that  arise  when  the  stress- strain 
curve  is  concave  upward.  These  comments  are  mainly 
due  to  only  one  of  the  authors  (HJH).  A  concave- upward 
curve  is  shown  in  Fig.  12.  The  smooth  curve  may  be 
approximated  by  two  or  more  chords,  as  shown,  joining, 
for  example,  the  points  0,  ( <r  i£i),  and  (  o*2,  £2) .  If  only  a 
single  chord,  joining  the  points  0  and  ( cr  2,  £2)  is  used,  a 
strain  wave  of  height  £2  will  travel  with  velocity 
The  equations  of  continuity  and  momentum  will  be  satis¬ 
fied,  and  the  equation  of  energy  will  also  be  satisfied  for 
the  linear  stress- strain  curve.  In  this  linear  case  it  is  well 
known  that  the  total  energy  is  half  kinetic  energy  and  half  strain 
energy.  But  strain  energy  is  found  by  integrating  under  the 
stress-strain  curve  from  0  to  £2.  If  we  use  the  linear  curve, 
energy  is  conserved  but  if  we  use  the  actual  curve,  or  the  two- 
chord  approximation  to  it,  we  find  that  the  potential  energy 
stored  in.  the  string  is  below  that  stored  in  the  linear  case.  If 
we  add  the  potential  and  kinetic  energies  we  find  that  they  do 
not  total  the  amount  of  work  done  in  generating  the  wave. 

A  way  out  of  this  dilemma  that  is  often  adopted 
is  to  assume  that  a  square  wave  of  height  £2  travels  with  a 
velocity  given  by  the  slope  of  the  single  chord,  and  to  assume 
that  the  work  done  in  generating  the  wave  goes  into  the  normal 
amount  of  kinetic  energy,  into  a  lesser  amount  of  strain  energy 
(corresponding  to  the  area  under  the  curve  that  is  concave 
upward)  ,  and  into  heat.  The  amount  of  heat  generated  is  just 
sufficient  to  conserve  energy,  but  the  way  in  which  it  is 
generated  is  usually  not  specified.  Some  insight  into  how  the 
extra  energy  might  be  converted  into  heat  is  given  by  two 
problems  solved  in  the  body  of  the  present  report:  the  problem 
of  the  equal  masses  joined  by  slack  strings,  and  the  problem  in 
which  a  disturbance  consisting  of  a  wave  of  fixed  height  was 
assumed,  with  an  oscillatory  wave  superposed  on  it.  In  each 
of  these  problems  it  is  possible  to  satisfy  the  equations  of 
momentum,  continuity,  and  energy. 
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A  possibility  exists,  however,  that  the  model 
usually  adopted  is  incorrect.  If  this  is  true,  the  treatments 
usually  given,  including  a  part  of  what  has  been  presented  in 
the  present  report,  do  not  correspond  to  reality.  They  are 
self-consistent  and  formally  correct  but  they  may  still  be 
physically  incorrect.  To  see  where  the  usual  assumptions 
may  be  wrong,  it  is  necessary  to  consider  the  nature  of  a 
wave  of  tension  as  it  propagates  into  a  string.  We  will  assume 
that  the  wave  is  a  square  wave.  Mathematically  the  wave  is 
infinitely  steep,  but  physically  it  must  certainly  have  a  finite 
thickness.  This  thickness  may  be  very  small,  but  it  is  still 
several  atomic  spacings.  Such  a  wave  is  shown  schematically 
in  Fig.  23  with  the  thickness  greatly  magnified.  Now  consider 
the  speed  of  propagation  of  this  wave  as  determined  from  the 
smooth  stress-strain  curve  of  Fig,  12. 

The  leading  edge  of  the  wave,  where  the  strain 
is  very  small,  must  propagate  with  the  velocity 

Ct  o  ~ 

where  the  subscript  0  indicates  the  origin  in  Fig.  12.  Where 
the  strain  is  higher,  the  propagation  velocity  is  higher  than  cq. 
Hence  by  the  usual  logic  a  wave  that  is  not  initially  steep  will 
become  steeper  because  the  regions  of  high  strain  travel 
faster  than  those  of  lower  strain.  However,  there  is  no  physical 
mechanism  that  will  permit  the  steepening  to  go  on  indefinitely. 

The  speed  of  propagation  of  a  small  strain  into  the  unstrained 
material  is  cq  and  the  entire  wave  can  propagate  no  faster. 

To  emphasize  this  point  we  will  mention  that  the 
elastic  properties  of  a  material  are  determined  by  the  forces 
between  the  various  atoms  and  molecules  that  make  it  up.  When 
the  material  is  unstrained  the  molecules  are  at  their  equilibrium 
distance  apart  and  a  small  strain  propagates  with  velocity  cq. 

Only  after  the  molecules  have  moved  fartier  apart  does  the  material 
have  a  higher  value  of  8  <r/  8  e  and  a  higher  strain  velocity.  In 
contrast,  if  we  were  dealing  with  a  compression  wave  rather  than 


(172) 
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a  tensile  wave,  the  speed  of  the  wave  would  not  be  directly 
determined  by  the  intermolecular  force  curves.  In  a 
compression  wave  there  is  no  practical  limit  on  the  magnitude 
of  the  force  that  can  be  applied  to  the  rear  surface  of  the  shock 
by  the  fully  strained  material.  The  pressure  simply  builds  up 
until  the  requisite  amount  of  force  is  applied.  Interatomic 
repulsive  forces  become  very  great  as  a  material  is  compressed 
and  so  there  is  no  practical  limit  to  the  process,  even  when  the 
particle  velocity  behind  the  wave  becomes  very  large.  In  a 
tensile  wave,  on  the  contrary,  a  force  that  is  too  large  will 
simply  pull  the  molecules  apart.  After  the  specimen  is  ruptured, 
no  further  work  can  be  done  on  the  unstrained  material. 

This  situation  is  reminiscent  of  what  happens  in 
a  gas,  where  a  compression  shock,  traveling  with  a  velocity 
greater  than  that  of  sound,  is  possible.  Analogous  expansion 
shocks  in  gases  do  not  exist.  In  expansion  we  find  only  expansion 
waves,  and  no  propagation  velocities  exceed  the  speed  of  sound. 
The  analogy  between  the  behavior  of  a  strain  wave  and  a  gaseous 
shock  wave  is  not  exact,  however,  and  should  be  used  with  care. 

If  the  arguments  given  above  are  correct,  no 
tensile  strain  wave  can  propagate  into  a  material  with  a  velocity 
greater  than  c^.  The  model  of  a  shock  wave  of  height  €z  travelirg 
with  velocity  c  =  Y cr  2/jpe z  as  required  by  eq  (51)  and  Fig.  12 
becomes  untenable  because  c>cQ.  Since  a  shock  wave  of  height  6z 
and  velocity  c  is  known  to  satisfy  both  continuity  and  momentum 
conditions  (but  not  energy)  do  we  not  make  matters  worse  by 
fixing  the  velocity  at  cQ  ?  We  do  not,  of  course,  if  the  forward 
propagation  velocity  actually  is  cq. 

The  problem  of  the  merging  of  two  finite  waves, 
one  overtaking  the  other,  was  treated  in  the  body  of  this  report. 
The  stress-strain  curve  for  that  problem  was  two  intersecting 
straight  lines  such  as  those  shown  in  Fig.  12,  and  the  wave 
heights  were  6i  and  Ez-€i  respectively,  so  that  each  pulse 
propagated  as  a  square  wave.  When  the  second  wave  overtook 
the  first,  the  propagation  velocity  of  the  merged  wave  abruptly 
changed  to  a  value  higher  than  that  of  the  leading  wave  and  lower 
than  that  of  the  overtaking  wave.  However,  on  the  basis  of  the 
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discussion  given  above,  we  expect  the  propagation  velocity 
of  the  merged  wave  to  be  the  same  as  that  of  the  leading  wave, 
with  a  sudden  change  in  strain  and  a  sudden  change  in  particle 
velocity  originating  at  the  point  of  overtaking.  In  order  to 
satisfy  the  conservation  equations,  a  reflected  wave  is  required. 
This  wave,  since  it  travels  in  strained  material,  will  have  the 
higher  propagation  velocity  of  the  second  wave. 

One  of  the  authors  (HJH)  made  several  attempts 
to  calculate  the  various  waves  and  show  how  continuity,  momentum, 
and  energy  conditions  could  all  be  satisfied.  He  did  not  get  a 
satisfactory  solution.  Whether  this  was  due  to  insufficient  mathe¬ 
matical  prowess  or  whether  the  problem  was  not  set  up  in  a  suit¬ 
able  form  is  not  known.  Most  of  the  assumed  solutions  that  were 
tried  out  began  by  assuming  waves  of  constant  amplitude  propagat¬ 
ing  in  both  directions.  It  may  well  be  that  waves  that  change  in 
amplitude  as  they  propagate  away  from  the  point  of  overtaking  are 
required.  Whatever  the  trouble,  there  seems  to  be  no  reason  why 
a  solution  of  the  overtaking  problem,  with  the  merged  wave  having 
a  propagation  velocity  of  cQ,  could  not  be  obtained.  This  should 
lead  to  an  under  standing  of  what  happens  when  the  stress- strain 
curve  is  a  smooth  curve,  concave  upward.  Then  comparison 
with  experiments  could  be  made. 

Some  attempts  were  made  to  solve  the  overtaking 
problem  involving  two  square  waves  by  computer  methods.  A 
stress-strain  curve  consisting  of  two  straight  lines  was  accepted, 
with  constants  chosen  so  that  the  second  wave  would  travel  twice 
as  fast  as  the  first.  Our  thought  was  that  we  could  plot  the  results 
and  determine  the  speeds  of  propagation  of  both  waves,  and  perhaps 
also  those  of  the  advancing  and  reflected  waves  formed  at  the  point 
of  overtaking.-'  The  attempts  were  unsuccessful.  Our  iterative 
method  was  able  to  deal  with  linear  stress-strain  curves  and  with 
curves  that  were  concave  downward.-  It  was  able  to  use  a  curve 
that  was  moderately  concave  upward  if  the  initial  strain  pulse  had  a 
gradual  rise,  presumably  because  the  concavity  steepened  the  prise 
without  making  it  entirely  square.  But  the  iterative  process,  when 
used  with  a  stress-strain  curve  that  was  concave  upward,  seemed 
always  to  fail  when  a  square  pulse  was  programmed  initially  and 
also  when  a  strain  pulse  that  rose  gradually  was  sharpened  by 
propagation. 
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